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1 Introduction 



Since the pioneering work of Bloch M, several questions concerning persistent currents 
have been answered. The conducting ring case has been largely discussed in the literature 
(2|, as well as the persistent current due to a point-like vortex Q. 

However, many questions remain open. For example, we might ask about standard 
relationships such as 

M = IV (2) 

where I, M, E, V, <J> are, respectively, the persistent current, magnetization, energy, area 
(denoted by V as volume in two dimensions) and magnetic flux through the system. 

These relations are clearly understood in the case of a ring. Do they apply to more 
general, for instance infinite, systems? More precisely, are they still correct, and if yes, in 
which case? 

The aim of this paper is to clarify these points in the case of two-dimensional systems. 
We will show that (|2|) holds true as soon as the system is invariant by translation and 
rotation. We will also insist on the role of the spin coupling, which happens to be crucial 
in the case of point vortices. We will examplify these considerations for an electron (or 
gas of electrons) coupled to 

- a point vortex superposed to an homogeneous magnetic field, 

- an uniform electric field and an homogeneous magnetic field (quantum Hall geome- 

try), 
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- a distribution of point-like vortices randomly dropped onto the plane according to 
Poisson's law, modelising disordered magnetic systems (analytical and numerical). 

The case of electrons on a ring threaded by a flux, and electrons in a plane coupled 
to a point vortex will be revisited in Appendix A. Some technical details on harmonic 
regularization for computing partition functions and the so called "Landau counting rule" 
will be given in Appendix B. 



2 Basic Definitions and Formalism 



Consider the two-dimensional quantum mechanical problem of an electron with Hamilto- 
nian (we set the electron mass and the Planck constant m e = h = 1) 

H= l -(p-eAf + V- e -a z B (3) 

V and A are, respectively, the scalar and vector potential (B = V x A) and a z /2 is the 
electron spin. The local current is (v = p — eA) 

-* 6 

j(f) = — (v\ r)(f\ + | f)(f\v) (4) 

The total magnetization is by definition 

M = \ J dr(r x <j(f)» • k + | (a z ) (5) 

where k is the unit vector perpendicular to the plane and ( ) means average over Boltz- 
mann or Fermi-Dirac distributions. In the Boltzmann case, one obtains the thermal mag- 
netization (Zp = Tr {e - ^}) 

Mp = ^-Tr{e-P H ((rxv)-k + a z )} (6) 
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We want to relate the thermal magnetization to a thermodynamical quantity. To this 
aim, we add a fictitious uniform magnetic field B' perpendicular to the plane. In the 
symmetric gauge, its vector potential is eA' = {e/2)B' x f. H becomes H(B') with the 
partition function Zp{B') = Tr |e _/3 ^( B ) j. First order perturbation theory in B' gives 

Zp(B') = Z (3 + e -^Tx {e-? H ((? x 0) • k + a z ) } + . . . (7) 
and thus, the magnetization, 

1 d\nZ 3 {B') 
Mp = — im 8 

M /3 -B'-»o dB' 

In the sequel, we will consider the spin-up and spin-down as two distinct physical situa- 
tions. It follows that in (|8|), the spin induced part of the magnetization is |cr 2 , thus the 
orbital part of the magnetization is 

«*. = i , im ?Msm . (9 , 

P p b'^o dB' 2 v; 

Note that in the particular case of H containing an homogeneous magnetic field B 
perpendicular to the plane, (||) narrows down to the usual formula 

IdlaZa 

= nn»r (10) 



and (^) to 

M orb_ I dlnZp e 

M ? -^~^B~~^° Z (A1) 



Note also that one can compute directly M°^ h using (||) by simply dropping the coupling 
of the fictitious B' field to the spin in Zp(B'), and, accordingly, using (0) in the case of 
H containing a homgeneous B field, the coupling of B to the spin in Zq. 

Let us now turn to the persistent current and consider in the plane a semi infinite line 
T> starting at r*o and making an angle o with the horizontal x-axis. The orbital persistent 



current I or (Fq,#o) through the line is 



(r — r*o) x (j(r) ) 

l orb (ro,0 ) = / d\r- r | — -A / . fc (1 2 ) 

jz> r — ro 



where — p^p^ • is the orthoradial component of j(f) to the line. It obviously depends 
on ro and 9q. Consider now systems rotationnally invariant around rb- / (Fq,#o) no 
longer depends on 9q and, without loss of generality, can be averaged over 9q. So 

I Orb (r ) = -L / dr ^-y^y^ • £ (13) 
2vr 7 |r — ro| z 

and for the Boltzmann distribution 



^)=^^i^"V^r' fc ) (14) 

Again, we want to establish a link between ( |l4|) and a thermodynamical quantity. We 
add to the system a fictitious vortex 0' located at ro the origin of the line, with potential 
vector and magnetic field 

eA v = a' "* (f :,; o) (15) 
eB v = 2ira'5(r — fo)k (16) 

where a' = <j>' /<j> -(fio is the flux quantum- can always be chosen in the interval [0, 1] by 
periodicity. The resulting Hamiltonian is 

H(a r ) = X -(v-eA- e A v f + V - ^a z (B + B v ) (17) 

with partition function Zp{a'). 

Though being a well-established fact ||], we recall here that the spin coupling term 
a z (B + B v ) in (^Jl7|) is necessary to define in a non ambiguous way point vortices quantum 
systems. The behavior of the wavefunctions at the location Fq of the fictitious point 



vortex (or of a physical point vortex contained in the original Hamiltonian (3)) has to 
be specified properly. This is precisely achieved in a scale invariant way with the spin 
coupling term, which materializes either a repulsive hard-core (spin-down) ira'5(r — fo) 
or attractive (spin-up) — ■na'5{f — fo) point contact interaction at fo. In the standard 
Aharonov-Bohm problem, the wavefunctions are chosen to vanish at fo, i.e. the spin-down 
repulsive prescription. However, we might as well take the attractive spin-up prescription. 
Still, the persistent current should clearly be insensitive to the nature of the fictitious 
vortex introduced to define it, and in particular on the way it is regularized at short 
distance. 

First order perturbation theory in a' gives 

Z p (a') = Z„ + a' (3 (tt je~^ (r Z^l^ ■ ^} + 7ca z Gp(r , f )) + . . . (18) 

where Gp(r,r') = (r\e~P H \r r ) is the thermal propagator for the Hamiltonian (||). We 
define the total current around fo as 

. . . e , dlnZgia') . , 

w^ssJa-arf <19) 

and, from (|T^), 

^-^^r^-h^ <20) 

Note that in the particular case of H containing a vortex <j> at fo, or consisting of a ring 
of center fo threaded by a magnetic flux (j) = afio, the orbital persistent current generated 
in the plane around fo, or in the ringF], (^) narrows down to 

b e dlnZp e Gp{fo,r ) 

^ {r » ) = 2^^r-r z Zp (21) 



If the first term in (|20| ) is usually discussed in the literature under the weaker form 
appearing in (pTI) valid for the point vortex and the ring P|-0|, the second term in (2C ,21) 



2 In the case of a ring, the propagator at the center of the ring is by definition nul. 
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triggered by the spin coupling to the propagator has been so far ignored. However, both 
terms are crucial, and if the spin term is absent, some contradictions do arise in the 
computation of orbital persistent currents. Indeed, orbital persistent currents -as well as 
orbital magnetizations- should vanish when H becomes free, and this is not the case if 
the spin term is ignored, as we will see later for the point vortex, or for the point vortex 
superposed to an homogeneous field. We will also encounter in the attractive point vortex 
case a situation where the orbital current by itself is not defined (in fact infinite), and 
where only the total (orbital + spin) current has a sensible physical meaning. 

To summarize, orbital magnetizations and currents are respectively given by @ and 
(pC|). They should both vanish when the original Hamiltonian H becomes free. Total 
magnetizations and currents are given by (Jsj) and (0). They obviously depend, as well as 
their orbital counterparts, on the spin degree of freedom defined in the original Hamiltonian 



If we consider systems that are both invariant by translation and rotation, we can 
proceed further by averaging (|T|) over -tq. Taking advantage of the identity || 



which generalizes the conducting ring situation to translation and rotation invariant two- 
dimensional systems. Furthermore, since for a translation invariant system Gp{rQ,fo) = 
Zp/V, (^) holds also for the total current and magnetization 




(!)• 




(22) 



we can write, using d6[p^). 




(23) 





(24) 
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Finally, (23.24) still hold for non rotational invariant systems if the persistent current I orb 
is understood as averaged over all directions 9q of the semi infinite straight line. 

To have access to the orbital magnetization and persistent currents of a gas of non- 
interacting electrons, Zg(B') or Zp(a!) have to be replaced in (p|,|2"0|), and Zg in ( JTT| , |2T| ) , 
by the corresponding grand-partition functions. For instance, (p~T| ) and ( pl| ) respectively 
become 

•"-JT-"?" (25) 

= ^lk^ M - l°^W-f»)f(H)} (26) 

where S / g(//) is the grand-partition function -ln!Hg(/i) = Tr |ln(l + e~ ^)|- for chemi- 
cal potential fx and f(E) is the Fermi-Dirac distribution. At zero temperature, the second 
quantized magnetization and persistent current at Fermi energy Ep are by definition, 
using the integral representation of the step function 9{Ep — H), 

^ roo giEpt 

M, F = -|^— Z,M,| W (27) 

Y roo giEpt' 

Ief = 2™ loo dt 'jr—; z ^y^> (28) 

where rj' , e' — > + . One should bear in mind that the number of electrons N may be fixed, 
implicitly determining Ep by 

^ roo e iE F t' 

N = ^iL dt '¥^? z ^ + < <29) 

We insist here that we regard the first quantized thermalized magnetizations and cur- 
rents (^[^), or (|8],[n^), as the basic objects from which their second quantized versions can 
be computed, for instance via (|27| ) and (^) at zero temperature. In particular if Mg or Ig 
happens not to depend on the temperature, that is to say if the magnetization or current 
for a quantum state do not depend on the energy of the quantum state, then the second 



quantized magnetization and current are simply given by M = NMp or I = Nip - see for 
instance in Appendix A the current for the point vortex case. 



3 Applications 



Z ? = —«— (30) 



3.1 The homogeneous magnetic field case 

As a warm-up exercise, consider an homogeneous magnetic field perpendicular to the 
plane. The spin-up (a z = +1) or down (a z = — 1) partition function 

V be a * b 
2tt(5 sinh b 

where b = (3lu c and u> c = eB/2 (eB > is assumed without loss of generality) leads to the 
orbital magnetization 

p P dB 2 2 \b J y ' 

which in this particular case does not depend on a z . 

Since the system is invariant by translation and rotation, we know a priori that the 
orbital persistent current is given by Ig rb = M^ rb /V . Let us however illustrate the con- 
siderations above by computing directly the orbital persistent current around a point 
r*o. Adding a vortex a' located at this point with a repulsive (spin-down) or attractive 
(spin-up) prescription, we get Q for < a' < 1/2 the partition function (see Appendix 
B) 

Z P {a) - Z P + 2^nhb{ a ~ 6 ~^nhb~ (32) 



One has also 



G£(?o,r ) = ±Z%* (33) 



Thus, ( pO] ) gives the orbital persistent current 

e dhxZZHa!) e e /l \ 

p 2vr/3a'^o da' 2V 2V \b J v ' 

which does not depend on r*o, and indeed satisfies Mp rh = I^ rb V. Note that it does not 

depend on a z , i.e. in particular on the short distance behavior of the fictitious vortex. 

3.2 The point vortex + a homogeneous magnetic field case 

We consider a point vortex carrying a flux 4>/4>o = a at the origin of the plane superposed 
to a homogeneous B field, and concentrate on the repulsive spin-down case, bearing in 
mind that the analysis in the attractive spin-up case would follow the same lines as in 
the point vortex case (see Appendix A). Using the partition function (^) with a' — > a, 



a z = —1, (11) gives the orbital magnetization in the spin-down case 



Mg rb = -(r-coth6)-^lfa-e-(^^ 



p - 2 l 6- CULnU; -^^V ^mhfe 

+-^-^- 1 > b (a - e^ bS ^)) +... (35) 
sinh b smh b J 

where . . . denote corrections which are subleading in volume. The vortex does not affect 

the leading volume term, which is therefore identical to the pure homogeneous magnetic 

field magnetization (]3l|). 

Let us now turn to the persistent current around the vortex location at the origin of 
the plane. 

Since, when a > 0, 

G^(0,0) = (36) 
1/3(0) = lp rb (0). From (21), the persistent current at leading order in volume reads 
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with indeed Mj? b / I^ b V . 

When a = on the other hand, the orbital persistent current should narrow down to 



the homogeneous magnetic field result (34). However, this does not happen for the current 



in pTj). But, if one pays attention to the fact that, when a = 0, 



^(0,0) = ^ (38) 

where Zn stands for the partition function of the homogeneous B field, then the spin term 
is such that in ( |2l| ) the orbital persistent current correctly narrows down to the orbital 
current for the homogeneous magnetic field, i.e. 



r-M\-T^*) + w-w<-i- a * h) (39) 



3.3 The case of the plane with an electric field 

Consider now a homogeneous magnetic field and an electric field E in the horizontal 
a>direction (Hall geometry) and ask how the magnetization and the persistent current 
([H|,34) are affected by the electric field. The spectrum is unbounded from below, making 



the use of thermodynamical quantities as the partition function hazardous. We propose 
to circumvent this difficulty by confining the electron in a harmonic well (for more details 
on harmonic regularization see Appendix B). The Hamiltonian in the symmetric gauge 
A = Bk x f/2 reads 

H = ^(p- eA) 2 - eEx + ^uj 2 f 2 (40) 
with indeed a bounded spectrum. Its propagator is 

G,(f,f') = Gf"\r- e-|,f' - e^)e^^^e^ (41) 
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where 

(42) 

is the propagator for the homogeneous B field in a harmonic well and u>t = \/w 2 + cJf. 
Note that in the thermodynamic limit uj —* 0, 

G«(r, r") = — ^ — exp ( — ^— — [cosh (3u> c (r — r') 2 + 2? sinh/?u; c (f x r") • fc 

2-7T sinh pw c L 2 sinh /xj c 



-sinh^ c (x + x')^ + ^£(-L - coth«^(^ + 2% - ,'))]} (43) 



still depends on f at coinciding points r = f ', a direct manifestation of the breaking of 
translation invariance. The partition function reads (without the spin coupling to the B 
field since we are interested in the orbital magnetization only) 

Zb = - - ei^) 2 (44) 

2 cosh (3ut — cosh j3ui c 

If E ^ 0, the partition function blows up in the thermodynamic limit uj —* 0. On the 

contrary, when E = 0, it correctly yields back the partition function for the homogeneous 

magnetic field if the proper thermodynamic limit prescription lim^^o (fjp" ^ ^ is taken. 

The orbital magnetization reads 

, e sinh 3uj c — ^ sinh Bait 
p 2 cosh Bujt — cosh Buj c 



which happens to coincide with (|3l| ) in the thermodynamic limit. The electric field has 
no effect on the magnetization of the system, a result that can be understood at the semi- 
classical level by the fact that the electric field only induces a translation of cyclotron 
orbits. 

Turning now to the persistent current around a point f*o, a fictitious vortex a 1 located 
at ro is added to the system, that we take to be repulsive. All what is needed is the 
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partition function of an homogeneous magnetic field in presence of an electric field and 
a harmonic well, and a repulsive vortex at rb, at leading order in a'. It can be obtained 
in perturbation theory along the lines developed in |?]]. The harmonic well has not only 
the virtue, as we have seen, to make tractable the computation of an otherwise diverging 
partition function with an E field, but it is also needed j^] to give an unambiguous meaning 
to the perturbative analysis in a', which is otherwise ill-defined. 

The first order perturbative correction to the partition function Zp(a') = Zp + a'Zp + 

Z} = 2/3 f dzdz^-^(d z - \uj c z)Gp{r, r%,^ (46) 

J Z Zq L 

where z = x + iy is the complex coordinate in the plane, gives an orbital persistent current 
equal to 

J^( r - ) = le-^k^^-^^JkVe^P^t / sinh/?Wc _ Wc sinhpuJ \ (48) 

M 2tt smhpuJt \ uJt J 

The current is exponentially damped when the point r*o around which it is estimated 
differs from the equilibrium position in the harmonic well r Q — e-jj? = 0. When E = 0, 
([Ig|) coincides, in the thermodynamic limit, with the orbital current for the homogeneous 
magnetic field (|34|), as it should. When E 0, this is still the case if tq = eE/uj 2 , 
which means that the effect of the electric field it to locate the current at the edge of 
the sample. However, as soon as tq ^ eE/ui 2 , the orbital curent Jg r6 (ro) vanishes in the 
thermodynamic limit. 
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4 The Magnetic Impurity Problem 



Consider || a random Poissonian distribution {f^, i = 1, 2, . . . , iV} of iV hard-core vortices 
at position r*j, carrying a flux (j) = a<p , where a can always be taken in the interval 
[0,1/2] without loss of generality. After averaging over disorder, the system is invariant 
by translation and rotation, so it is sufficient to compute M orb to get the orbital persistent 
current, since they are related by (§). The Hamiltonian is 

H = ~(p-eA(r)) 2 + ^B(r) (49) 

with 

f^W = (50) 

N 

eB(r) = 2-Ka^5(r-ri) (51) 



=i 



B(f) might be replaced by its mean value 



(B) = *m (52) 

e 

a mean field approximation that is valid in the limit a — > 0, and p = N/V, the vortex 
density, fixed ||. 

In a Brownian motion approach ||, for a given configuration of vortices, the partition 
function is 

Zp = Z° p ^exp ^ Znriicxj ^ (53) 

where ( ){c} means averaging over all closed Brownian curves of length (3, and rii is the 
winding number around the vortex at location fj. Z% = V/(2nf3) is the free partition 
function. 
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Averaging over all vortices configurations leads to 

Zp = Z° p ^exp \p-£ S n (e 2 *™ n - 1)) J (54) 

where S n stands for the arithmetic area of the n-winding sector of the given Brownian 
curve chosen in {C}. Due to scaling properties of Brownian curves ^, the random 
variable S n scales like /?, thus the rescaled variables 

2 1 

S = — ^2 sin (nan) A = — ^ S n sin(27ran) (55) 

P n P n 

The probability distributions of S and A are actually independent of (3 with (S)r C ~> = 
na(l — a) and (A)r C \ = 0. 

To get the magnetization, we add a fictitious magnetic field B' , with partition function 
(without the spin coupling to the fictitious B' field since we are interested in the orbital 
magnetization only) 

Zp(B') = Z° p L-P P {S-iA)+iPB>A\ (56) 

where A = (1//3) J2n n ^n is the rescaled algebraic area enclosed by the Brownian curve in 
{C}. From (||), and since S — > S, A — > —A, and A — > —A when one spans a Brownian 
curve in the opposite direction, the orbital magnetization is 

(e-PP s sm(f3pA)A) snx 

T orb _ \ {JO} 

(e'^cos((3 P A)} {c} 



M T = ~ ,„-aoS^ ta „^ X t ( 57 ) 



Mp rb is actually only a function of (3p and a, odd in a. Thus, for a 6 [0, 1/2], necessarily 

oo 

M% rb = (1 - 2a)F(f3p, a(l - a)) = (1 - 2a) £ (flo)" £ a mn (a(l - a)) m (58) 

n=l m>n 

which can in principle be obtained in perturbation theory || ||]. 
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The p n a m term describes an electron interacting m times with n vortices, so m > n. 
The mean magnetic field ( |52|) scales like pa, therefore the mean field terms are those with 
m = n. They are in fact easy to obtain, since a — > corresponds to the mean field regime. 



Thus, from (31,52 



(M)- 6 = |(^-coth{6)) (59) 

where (b) = f3(u> c ), with (lu c ) = e(B) /2 = irpa. From (|58|), one infers that M| rb neces- 
sarily contains 

Mf\ mean = (1 - 2a) | {^j - coth {b)'j (60) 

where (b)' = (3irpa{\ — a). 

It is even possible to go a little bit further, by observing that the p n a n+1 terms come 
from the homogeneous (B) field + one vortex case ||. The first correction to the mean 
field approximation (BO) is 



P =M P \mean+ea(l-a)(l-2a) + - _ e _ m , + (1 _ e _ 2(fe) , )2 J+--- 

(61) 



In Fig. 1 {(3p = 1), the agreement is rather good between (61) - full curve - and the 
numerical simulations - points - based on (|57|). We generated 2000 random walks of 100000 
steps each one. However, the situation becomes less transparent for higher (3p values. 
Clearly, the perturbative analytical approach need more and more corrections coming 
from (B) + two vortices, .... The numerical simulations indicate that, for < a < 1/2, 
Mp rb is negative whatever (3p is. Even if it is physically clear why it should be so since 
the vortices enforce a clockwise rotation of the electrons when < a < 1/2, an analytical 
proof is missing. 
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A Appendix A: Persistent current for the ring and the point 
vortex 



A.l The case of the ring revisited 

Consider a ring of radius R threaded by a homogeneous magnetic field B perpendicular 
to the ring, with flux through the ring. Let us denote as usual a = ec/>/27r = <j)/4> 
the value of <j> in unit of the quantum of flux. The system is periodic in a with period 
1. Reversing the sign of B does not change the physics either -the partition function is 
symmetric around a = 1/2, whereas the magnetization and current are antisymmytric-, 
so one can always restrict < a < 1/2. Clearly, the orbital persistent current changes its 
sign with <j). But (j) = ±1/2 are identical, thus necessarily I°p h (a = 1/2) = 0. 

To summarize, both first and second quantized orbital currents should be defined in 
the interval a S [0, 1/2], and due to obvious symmetry considerations, they should vanish 
when a = 0, 1/2. 

In unit 2R 2 = 1 the partition function is 



m=oo 

\2 



Z p {a)= Y, (62) 



m=— oo 



The thermal current is formally defined as in (^). However, it is equivalent to thread 
the ring with an homogeneous B field or at its center with a vortex line of flux (j). It follows 
that (p0|) narrows down to (|2l|). Clearly no distinction has to be made either between the 
orbital and total magnetization (or current) since by geometry the propagator at the center 
of the ring is zero. Thus 

r _ rorb _ e dhlZ^ 

1/3 ~^p^r (63) 

and 
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Since cj) = a2Tr/e = BnR 2 , one necessarily has Mg = IgirR 2 . 



One obtains 



I 



sinh(2/c/3a) k 



sinh(fe/3) 

or equivalently, using the reciprocity formula, 



Zg(a) 



(65) 



(66) 



-y 



sin(27rA;a) 



-1) 



Clearly, from (65,67) 



■' ^ sinh^vrV/?) 
I 3 (a = 0) = = Ip{a = 1/2) 



(67) 



(68) 



as it should. Note that (65.67) are, to our knowledge, the simplest expressions obtained 
so far for the thermalized persistent current on a ring. 

© is the sine Fourier series of a 
periodic and odd function, making the symmetries of lg explicit. 

Considering now a second quantized gas of electrons on the ring, starting from Ig, the 
second quantized persistent current at zero temperature and Fermi energy Ep is given by 
(p£|). However, since individual persistent currents for each quantum states are somple to 
obtain, one may sum individual currents up to the Fermi energy^. The current carried by 
the state \m > of orbital momentum m and energy (m — a) 2 is I m = §(w — a). It follows 
that if N is even, 

Ie f = ^N(l - 2a) 



2vr 



and if N is odd, 



L E F 



-~aN 

7T 



(69) 
(70) 



From now on we revisite the approach followed in 
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These results seem misleading since the currents in ( p9|j70| ) do not vanish either at a = 
or a = 1/2. However, their derivation in terms of individual quantum states is valid for 
a spectrum with no degenerate states, i.e. if and only if < a < 1/2. Indeed, when 
a = 0, the states |m > and | — m > are degenerate, whereas when a = 1/2, the states 
|1 — m > and \m > are degenerate. It follows that ( |69| ) is uncorrect when a = 0, and ( |70| ) 
is uncorrect when a = 1/2. If the degeneracy is properly taken into account, meaning that 
a given degenerate doublet of states contributes to the current as Jm + J - m when a = 0, or 
h- m +i m w ken a _ \/2, the total current ends up to vanish as it should when a = 0, 1/2. 

We conclude that the persistent current for a gas of electrons at zero temperature on 
a ring threaded by a flux a<p is respectively given by ( |69D or ( |70D when ./V is even or odd 
if < a < 1/2, and if a = 0, 1/2. 

A. 2 The case of point vortex revisited 

Consider finally a point vortex of flux <ft = a<pQ at the origin of the plane, which can 
be either hard-core (spin-down repulsive prescription), or attractive (spin-up attractive 
prescription), clearly two different physical situations. 

As in the magnetic impurity and ring cases, one can always restrict the interval of 
definition of a £ [0,1/2]. However, to illustrate the non trivial role of the spin, we will 
explicit the results in the whole interval a £ [0, 1[. 

The spin-down and (up) partition function respectively read (see Appendix B) 

Z~ - Z° p = ^a(a - 1) < a < 1 (71) 

Z+-Z$ = \a{a + l) 0<a<^ 

Z+-Z% = l( a -l)(a-2) \< a < 1 (72) 
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leading to the total current around the originF] 



Ip = y (« - I) < a < 1 (73) 



It = \<a<l (74) 

which is discontinuous at a = 1/2 in the spin-up case. 

We would like to identify the orbital part if of the currents (73,[fj|). 
If a = 0, i.e. the free case, by definition 

G+(0,0) = G£(0,0) = ^ (75) 



One rightly deduces 

e . 1, e 

_ — (a + <7 Z — — &z — 



if = lim ^(q + ,4) - a, ^ = (76) 



as it should. 

If < a < 1, in the repulsive spin-down case, 

G^(0,0) = (77) 

We conclude that the orbital persistent current generated by a repulsive vortex at the 
origin is 

if- =0 a = (78) 



if' = ^{a-\) 0<n < 1 (79) 



4 Computing the current around a point fo would require the partition function of a vortex a at the 
origin and another fictitious vortex a 1 at f*o at leading order in a' . 
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On the other hand, in the attractive spin-up case, the propagator at the location of 
the vortex is infinite as soon as a / 0, thus I^ rb+ diverges. We conclude that the orbital 
persistent current generated by an attractive vortex is not defined (in fact infinite), only 
the total current (orbital + spin) is defined and given by (f74|). 

The thermalized persistent currents ( ff^JT^ ) do not depend on the temperature. It 
follows trivially that for a gas of N electrons, in the interval < a < 1/2 for instance, the 
persistent current rewrites as 

T» =e Pe (a + ^) (80) 

where p e = N/V is the electron density. 

The magnetization induced by a vortex in the plane (|l0| ) is deduced from the partition 
function of a vortex a + a fictitious B' field, i.e. (|32]), with a' — > a and B — > B' . One gets 

M or b - = _^^l a[l _ a){ }_ _ a) Q<a<1 (81) 

and 

M „ rb+ = _^ a(1 + a)( l +a) <«<i 

M ^ b+ = _^ (1 _ a)(2 _ a)( 3_ a) i<a<l (82) 

which are subleading in volume compared to the homogeneous field case. The magnetiza- 
tion for a gas of electrons is, for < a < 1/2, 

M orb,a 2 = _£ a((J + a )(— + a)— i (83) 

6 v A 2 'e-Pi* + 1 

with the additional Fermi factor at zero energy If it is to be compared to the current I Uz 
in (|^), one sees that both the current and the magnetization induced by a point vortex 
in the plane for a gas of electrons have the same leading volume behavior. 
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B Appendix B: Harmonic regularization and Naive Landau 
counting rule 



B.l Harmonic regularization 

Partition functions for systems with, in the thermodynamic limit, a continuous spectrum 
(point vortex) or a discrete but infinitely degenerated spectrum (homogeneous magnetic 
field) can be obtained in a non ambiguous way by regularizing the system at long distance 
by adding an harmonic term uP'r 1 /2 to the Hamiltonian, computing the partition function 
in the presence of the regulator, and then taking the thermodynamic limit uj — » with 
the prescription l/(/?u;) 2 — ► V/X 2 , where A = \/2Wj3 is the thermal wavelength. 

Let us illustrate this point by computing the partition function for the repulsive point 
vortex (|7l|). In the absence of a harmonic well, the spectrum is the continuous free 
spectrum with for eigenstates Bessel functions with shifted angular momentum m — > 
m — a. It is quite tedious and subtle to compute the partition function directly in the 
thermodynamic limit S. Consider rather the same problem in an harmonic well. The 
spectrum is now discrete with a finite degeneracy 



with n > 0. This is nothing but the spectrum of a two-dimensional harmonic oscillator 
with shifted angular momentum m — > m — a. The partition function reads 



E, 



'n,m — 



u>(2n + I'm — a\ + 1) 



(84) 



cosh/Ma-1/2) 
P 2sinh/3u;sinh/?u;/2 

Of course it diverges when uj — > 0, but the difference Q 



(85) 



ryiO ryOLJ 

A /3 ~ A /3 — 



sinh/3a;a/2 sinh/3a;(a — l)/2 
sinh (3ui sinh (3u> / 2 



a(a — 1) 



2 



(86) 
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remains finite^] in the thermodynamic limit and yields (fH|). 

The partition function for an homogeneous magnetic field ( |30| ) can be obtained along 
the same lines, with the spectrum 

E nm = u! t {2n + \m\ + 1) - oj c m (87) 

and in the case repulsive point vortex + homogeneous B field, with the spectrum 

E nm = uj t (2n + \m — a\ + f ) — iv c (m — a) (88) 

In the latter case, one computes Z^(B,a) — Zg(B,0) to obtain, in the thermodynamic 
limit u -»• 0, (H). 

B.2 Naive Landau counting rule 

The so-called "Landau counting rule" in terms of individual states happens to lead to 
sometimes uncorrect results, although one sums only over a finite number of Landau 
quantum states. Again, correct results are recovered if the harmonic regularization is used 
instead. 

Consider for example the spectrum (|88|) where, in the thermodynamic limit oj — > 0, i.e. 

cot — > lo c , the vortex induces a modification of the Landau spectrum only for a finite number 

of states of energy tv c (2j + 1 + 2a), j > 0, namely the states of orbital angular momentum 

—j<m<0. Computing the partition function ( |88|) directly in the thermodynamic 

limit with this shifted Landau spectrum would lead to a result different from (p2|). The 

harmonic regulator adds a contribution to the partition function from the states m > in 

(|SS|), that remains finite in the thermodynamic limit, eventhough it is an infinite sum of 

vanishing terms, and yields the correct expression (|32[). 

The same result can be obtained by Zeta function regularization - see the second reference CMO in 

I- 
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One can illustrate the failure of the "Landau counting rule" on an other example. 
Consider the thermalized orbital persistent current for an homogneous magnetic field (|34|) 



-, oo oo 
n 

M n=0 m=— oo 

in terms of the individual persistent current I nm for the Landau quantum states | n, m ) 
of orbital momentum m and energy -E njm = (2n + 1 + |m|)o; c — mu> c (i.e. (^7j) in the 
thermodynamic limit). The individual currents write 

rorb _ e / i u 1 v ew c m 

7 «rm= ^(n,m\-\n,m) = —(—-1) if m + 

= —(-1) ifm = (90) 

Despite the fact that the number of states in a given Landau level is infinite, again only 
a finite number of states contribute to the persistent current for each Landau level. Thus 
the sum ( |SS| ) is convergent. Yet, following this procedure, one gets the uncorrect result 
Ip rb = — 2yCoth6. This is again due to the fact that, in this naive "Landau counting 
rule", if the m > states yield vanishing individual currents, they are still an infinite 
number of them. To give an unambiguous meaning to this infinite summation of vanishing 
contributions, a long distance regulator is needed, bearing in mind that the thermodynamic 
limit should be taken afterwards the summation over the quantum numbers n, m. Indeed, 
if one adds the harmonic regulator to the Landau Hamiltonian (the spectrum is now given 
by (|87l)), one finds 



2tt \m\ <jJt 
ecu, 



(-1) ifm = (91) 

where ujt = %/uS 2, + uS*. Now the m > do have individual non vanishing persistent 
currents, and their total contribution to the thermalized current does not vanish in the 
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thermodynamic limit. It is precisely equal to 2y^, a contribution which, when added to 
those of the m < states, yield the correct result ([34|) for Ig rb . 

The same situation happens for the thermalized current induced in the plane by a 
repulsive point vortex superposed to an homogeneous B field with individual currents 

Z7r m — a uv 



Acknowledgements 

We would like to acknowledge useful conversations with E. Akkermans, G. Montambeaux, 
R. Narevich and D. Ullmo. 



References 

[1] F. Bloch, Phys. Rev. 137 (1965) A787 

[2] L. P. Levy, G. Dolan, J. Dunsmuir and H. Bouchiat, Phys. Rev. Lett. 64 (1990) 
2074; H. F. Cheung, Y. Gefen, E. K. Riedel and W. H. Shih, Phys. Rev. B 37 (1988) 
6050; H. F. Cheung, E. K. Riedel and Y. Gefen, Phys. Rev. Lett. 62 (1989) 587; M. 
Buttiker, Physica Scripta T54 (1994) 104; M. Buttiker, Y. Imry and R. Landauer, 
Phys. Lett. 96A (1983) 365; Y. Avishai and M. Kohmoto, Ben Gurion University 
Report for a review and references see R. Narevich, Technion Haifa Thesis 5757 (july 
1997) 

[3] E. Akkermans, A. Auerbach, J. E. Avron and B. Shapiro, Phys. Rev. Lett. 66 (1991) 
76; A. Comtet, A. Moroz and S. Ouvry, Phys. Rev. Lett. 74 (1995) 828; A. Moroz, 
Phys. Rev. A53 (1996) 669; for a review and references see R. Narevich, Technion 



25 



Haifa Thesis 5757 (july 1997); see also Yu. A. Sitenko and A. Yu. Babansky, Bo- 
golyubov Institute Report (1997) 

[4] A. Bergmann and G. Lozano, MIT Report (1993); A. Comtet, S. Mashkevich and S. 
Ouvry, Phys. Rev. D52 (1995) 2594 

[5] J. Desbois, C. Furtlehner and S. Ouvry, Nucl. Phys. B453 [FS] (1995) 759; J. Phys. 
I France 6 (1996) 641 

[6] J. Desbois, S. Ouvry and C. Texier, Nucl. Phys. B500 [FS] (1997) 486 

[7] A. Dasnieres de Veigy and S. Ouvry, Nucl. Phys. B388 [FS] (1992) 388 

[8] A. Comtet, J. Desbois and S. Ouvry, J. Phys. A : Math. Gen. 23 (1990) 3563; W. 
Werner, These, Universite Paris 7 (1993) 

[9] A. Comtet, Y. Georgelin and S. Ouvry, J. Phys. A : Math. Gen. 22 (1989) 3917 



26 




-0.05 1 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 



Figure 1: The orbital magnetization in unit e = 1 in the magnetic impurity problem for 
(3p = 1. Comparison between the analytical computation (|6l|), full curve, and numerical 



simulations (57). The points were obtained by generating 2000 curves of 100000 steps 
each one. 



27 



